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Abstract. We present recent results on charmonium and charmonium-like states from the BaBar 
B-factory located at the PEP-II asymmetric energy e + e~ storage ring at the SLAC National Accel- 
erator Laboratory. 
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INTRODUCTION 

The charmonium spectroscopy has been revitalized by the discovery of many new 
states above the open charm threshold. While some of them appear to be consistent 
with conventional charmonium others don't seem to behave like standard meson and 
could be made of a larger number of constituent quarks. While this possibility has been 
considered since the beginning of the standard model, the clear identification of such 
states would represent a major revolution in the picture of the quark model. The possible 
interpretations of these states beyond the conventional mesons are: 

• molecules: these are loosely bound states of a pair of mesons ([Qq] [q Q\) [1]. This 
system would be stable if the binding energy were to set the mass of the states 
below the sum of the two meson masses. The dominant binding mechanism of this 
state is pion exchange. 

• tetraquarks: a bound state of four quarks usually represented as ([<2<?][</<2]) [2]. 
Strong decays proceed via rearrangement processes. This interpretation predicts 
many new states and the existence of states with non-zero charge. 

• hybrids: a bound state of a quark- antiquark pair and a number of constituent gluons. 
The lowest-lying state from lattice calculation is expected have mass about 4.1 
GeV/c 2 [3]. While states with exotic quantum number such (0 + ~, 1~ + , 2 + ~) 
are unambiguos hybrids other states can be distinguished from the conventional 
hadrons by characteristics of their decay process. 

In addition there is the possibility that some observed enhancements are simply due to 
threshold effect: a given amplitude might be magnified when new hadronic final states 
become energetically available. 

Even if several experiments [4] have performed many measurements of these 
charmonium-like states the overall picture is not clear. In the next section we will 
review recent BABAR results on this field. 



CHARMONIUM PRODUCTION @ THE B -FACTORIES 



The B -factories are an ideal place to study charmonium since charmonium states are 
copiously produced in a variety of processes: 

• B decays: B mesons decay to charmonium in about 3% of cases. Charmonium 
states of any quantum number can be produced. 

• Double-charmonium production: in this process, observed for the first time by 
Belle [5], a J /y or a y(2S) is produced together and exclusively with another 
charmonium state. 

• 77 fusion: two virtual photons are emitted by the colliding e + e~ pair (e + e~ — > 
e + e~Y*Y* ~~ ^ e + e~(cc)). States with C=+l are formed. 

• Initial state radiation (ISR): where a photon is emitted by the incoming electron or 
positron (e + e~ — > y(cc)), only states with J PC = 1 are formed. 

STUDY OF THE J/\}/(Q IN TWO-PHOTON INTERACTIONS. 

The Y(3940) was observed for the first time by Belle [6] in B decays and then confirmed 
by BABAR [7]. In a re-analysis [8] of the BaBAR data sample the precision of the Y(3940) 
parameters was improved and in addition evidence for the decay X(3872) — > J /yea 
was found. This confirmed an earlier Belle claim [9] for the existence of this decay 
mode. A subsequent Belle paper [10] reports the evidence of a structure in the process 
77 — > J/yrco that they dubbed the X(3915) with mass and width values similar to those 
obtained for the Y(3940) by BABAR [7]. BAR4Rhas recently performed a study of the 
process YY^-J/yrco [11] to search for the X(3915) and the X(3872) using a data sample 
of 519 fb _1 . Figure 1 shows the reconstructed J/ya) mass distribution after all the 
selection criteria have been applied. A large peak at near 3915 MeV /c 2 is observed 
with a significance of 7.6 <7. The measured resonance's parameters obtained from a 
maximum likelihood fit are m x(3915 ) = (39 1 9.4 ±2.2± 1.6)MeV/c 2 , r x(3915 ) = (13 ± 
6 ± 3) MeV. The measured value of the two-photon width times the branching fraction, 
is r yr (X(3915) x <^(X(3915) ->■ J/\\f(o) = 52 ± 10 ±3 eV and (10.5 ± 1.9 ±0.6) eV 
for the spin hypothesis J=0 and J=2, respectively, where the first error is statistical and 
the second is systematic. In addition a Bayesian upper limit (UL) at 90 % confidence 
level (CL) is obtained for the X(3872): r yr (Z(3872) x ^(X(3872) -> J /yd) < 1.7 eV, 
assuming J=2. 

STUDY OF THE T] c 7r + 7r IN TWO-PHOTON INTERACTIONS. 

Studies of charmonium like-states in recent years have been performed using mainly the 
J I \irn + n~ final state, no search using the r\ c K + K~ final state has been conducted. Such 
a search may shed light on the quantum numbers of the charmonium-like-particle or the 
internal dynamics of these states. In particular, it has been suggested that if the X(3872) 
is the \ x Di state r/ C 2 then the branching fraction ^(X(3872) — > r\ c K + K~) could be sig- 
nificantly larger than ^(X(3872) ->■ J/yK + K~). The quantum numbers J PC = 2~+ of 
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FIGURE 1. The efficiency-corrected invariant mass distribution for the J/y/co final state. The solid line 
represents the total fit function. The dashed line is the background contribution. The solid histogram is 
the non J/y/co background estimated from sidebands. The vertical dashed line is placed at the X(3872) 
nominal mass. 
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FIGURE 2. Distribution of (a,c,e) m(K®K + n n + n ) with the fit function overlaid for the fit regions 
of the (a,b)% c2 (\P), (c,d)Vc(2S), and (e,f) X(3872), X(3915) and Xci{2P). The vertical dashed line in (f) 
indicate the peak mass positions of the X(3872), X(3915), and Xdi^P)- 



the r\ C 2 are consistent with the results of an angular analysis of X(3872) — > J /i//7r + 7T~ 
[12] and would allow production of X(3872) in two-photon fusion. BaBAR using 474 
fb _1 studied the process JJ— YX — > r} c n + n~[13] where X stands for one of the reso- 
nance Xc2^P), r) cl (2S), X(3872), Z(3915) or Xdi^P). The r\ c was reconstructed via 



its decay to K®K + n , with K® — >• k + k . The signal yield for each X resonance is ex- 
tracted from a two-dimensional fit to m(K®K + n~) and m(KsK + n~n + n~). Figure 2 
presents the two dimensional fits around each of the resonances. No significant signal 
is observed in any of the fits. Table 1 summarizes all the fit result. Upper limits are 
obtained on the branching fraction SS{r\ c {lS) ->■ t) c k + k~) < 7.4% @90%C.L.) and 
&(Xd2(lP) ->• r] c K + K-) < 2.2% @90%C.L.) 

TABLE 1. Results of the yy ^ rj c n + n^ fits. For each reso- 
nance X, we show the peak mass and width used in the fit and 
the 90% CL upper limit on the product of the two-photon par- 
tial width Tyy and the X — > T] c %% branching fraction. For the 
X(3872) and the X(3915) we assume J = 2. 



Resonances 


M x (MeV/c 2 ) 


T x (MeV) 


UL (eV) 


Xcl(lP) 


3556.20 ±0.99 


1.97±0.11 


15.7 


r? c (25) 


3638.5 ±1.7 


13.4±5.6 


133 


X(3872) 


3871. 57 ±0.25 


3.0±2.1 


11.1 


X(3915) 


3915.0±3.6 


17.0±10.4 


16 


Xa{2P) 


3927.2 ±2.6 


24±6 


19 



SEARCH FOR THE Z\ (4050)+ AND Z 2 (4250)+ 

In 2008 the Belle Collaboration reported the observation of a resonance-like structure 
called the Z(4430) + decaying to y(2S)n + in the study of the process B ->■ y(2S)Kn 
[14]. This claim generated a great deal of interest [15], [16] since such states must have 
a minimum quark content ccdu, and thus would represent an unequivocal manifesta- 
tion of four-quark meson state. The BaBAR collaboration searched the Z(4430) + in a 
similar analysis of the Belle collaboration in the process B — > y/(2S)Kn and also in 
B — > J/xyKn [17] but they did not find any structure neither in the J/y/K nor in the 
\j/(2S)n. Recently Belle performed an amplitude analysis of the J/y/Kn [18] system 
finding no significant evidence of the Z(4430) + in agreement with BABAR . In 2009 the 
Belle Collaboration reported also the observation of two resonance-like structures sim- 
ilar to the Z(4430)+ in the study of the 5° ->■ XaK^K + [19]. These new structures 
were labeled as the Zi(4050) + and Z 2 (4250)+, both decaying to XdX + . BaBAR using 
a data sample of 429 fb _1 has recent searched for the Zi(4050) + and Z2(4250) + in the 
process B° ->■ Xc\K~n + and in the decay B + ->■ K®Xcin + [20]; with the Xci -> J/WY- 
In the BABAR analysis the Xc\X + mass distribution, after background subtraction and 
efficiency-correction, has been modeled using the angular information from the Kit 
mass distribution that has been represented using the Legendre polynomial moments. 
The excellent description of the Xd^ + mass distribution given by this analysis ap- 
proach shows that there is no need for any additional resonance to model the distri- 
bution. Figure 3 shows the result of the fit done on the Xd^ + mass spectrum using 
two or one scalar Breit-Wigners with parameters fixed to the Belle measurement. In 
all the fit cases there are no significant resonant structure, the statistical significance 
obtained is very low < 2c. The ULs on the 90 % CL on the branching fractions are 
for the one resonance fit: ^(5° -> Z + K~) x ^(Z+ ->■ Xc\^) < 4.7 x 10~ 5 while 
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FIGURE 3. (a),(b) Background-subtracted and efficiency-corrected %c\% mass distribution for B — > 
%c\Kn. (a) Fit with the Z\ (4050)+ and Z 2 (4250) + resonances, (b) Fit with only the Z\ (4050)+ resonance. 
(c),(d) Efficiency-corrected and background-subtracted %e\% mass distribution in the K% mass region 
where Belle found the maximum resonance activity: 1.0 < m 2 (Kn) < 1.75GeV 2 /c 4 . (c) Fit with 



Zi (4050)+ and Z 2 (4250)+ resonances, (d) Fit with only the Z(4150)^ 
indicate the fitted resonant contributions. 



resonance. The dot-dashed curves 



for the two resonances fit are: 3§(B° -)> Z\K~) x &{Z\ -)> % C \TZ + ) < 1.8 x 10~ 5 and 
B8(B° Z+K-) x ^(Z+ -»■ ^ c i7T + ) < 4.0 x 1(T 5 . 

STUDY OF THE J/yn+ir SYSTEM VIA INITIAL STATE 

RADIATION (ISR). 

In 2005 BaBAR discovered the Y(4260) in the process e + e~ ->■ (4260), with the 
y(4260) — > J /\yiz + Tl~ [21]. Since it is produced directly in e + e~ annihilation it has 
J PC = 1 . The observation of the decay mode J/yn :°7T [22] established that it has 
zero isospin. However it is not observed to decay to D*D* [23], nor to D*D* [24], 
so that its properties do not lend themselves to a simple charmonium interpretation, 
and its nature is still unclear. A subsequent Belle analysis [25] of the same final state 
suggested also the existence of an additional resonance around 4.1 GeV/c 2 that they 
dubbed the Y(4008). BAR4Rhas performed recently a new analysis [26] of this pro- 
cess using 454 fb~ l . In this new study the region below 4.0 GeV / c 2 has been studied 
for the first time. As shown on Figure 4(a) in that region there is an excess of events 
above the sidebands background. To understand the nature of this contribution a 
detailed study of the ty/(2S) line shape has been performed and the results is that it 
is not possible to discount the possibility that its due to J/\yn + n~ continuum cross 
section in this region. Figure 4(a) shows the fit to the mass distribution. A 

clear signal for the Y(4260) is seen; the values obtained from an unbinned-maximum- 
likelihood fit are: m r(4260) = 4244 ±5 ±4 MeV/c 2 , r y(4260) = 114+J^ ±7 MeV and 
T ee x ^(J/\j/7t + 7t-) = 9.2 ± 0.8 ± 0.7eV. There is no evidence for the Y(4008) seen by 




FIGURE 4. (a): the J jx\fK + K~ mass spectrum from 3.74 GeV /c 2 to 5.5 GeV /c 2 ; the point represent the 
data and the shaded histogram is the background from the J/y/ sideband; the solid curve represent the fit 
result, (b) the n + n~ distribution of the from the Y(4260) decay to J/\j/n + n^ . The solid curve represent 
the fit using the model described in the text. 
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FIGURE 5. The y/(2S)7T + 7r~ invariant mass distribution from threshold to 5.95 GeV jc 2 for the 
V(25) — > //y/7r + 7r~; the points with error bars represent the data in the y(2S) signal region, and the 
shaded histogram is the background estimated from the yf(2S) sideband regions. The solid curve show the 
result of the fit. 



Belle [25]. 

In the new BaBAR analysis a detailed study of the system from the Y(4260) de- 

cay to J /\yiz + K~ has been performed. The 7r + 7T~ mass distribution shown on Figure 
4(b) seems to peak around the /o(980) mass; however the peak is displaced from the 
nominal /o(980) position, since it is around 940 MeV /c 2 . The fact that the peak is 
displaced and the particular shape of m(n + n~) distribution seems to suggest a pos- 
sible interference between the /o(980) and m(n + iz~) continuum. To test this possi- 
bility the /o(980) line shape is taken from the BaBAR analysis of the D$ — > k + tz~k + 
[27] and this amplitude as been used in a very simple model to describe the n + 7l~ 
mass distribution: | y/ pol + e'^F^g^ \ 2 where pol is a polynomial function used to de- 
scribe the m(n + Ti~) continuum and ^(980) is the amplitude from the D$ — > iz + k~k + 
[27] analysis; <j) allows for a phase difference between these amplitudes. The result 
of this study is shown on Figure 4(b) and indicates if there is a real /o(980) con- 
tribution to the decay of the Y(4260) to J/\yn + n~ its contribution is not dominant: 
^(y4260->^/y/Q(980),/o(980)->g+a-) _ n? , 



STUDY OF THE \j/(2S)7i + 7i- SYSTEM VIA INITIAL STATE 

RADIATION (ISR). 

In addition to the Y(4260), two more J PC = 1 states, the Y(4360) and the Y(4660) 
have been reported in ISR production e + e~ — > \y(2S)n + n~ [28], [29]. While the 
Y(4360) was discovered by BABAR [28] and then confirmed by Belle [29] the Y(4660) 
was only observed by the Belle Collaboration. BABAR performed a new analysis using 
all its available dataset collected ad the Y(nS), n=2,3,4; that corresponds to an integrated 
luminosity of 520fZ? -1 . The \i/(2S)n + n~ mass spectrum for the y(2S) ->■ J/yK + K~ 
is reported in Figure 5 BABAR observes two resonant structures, that have been in- 
terpreted as the Y(4360) and the Y(4660), respectively. The parameters values ob- 
tained from an unbinned-maximum- likelihood for the first resonance are m Y (^6Q) = 
4340 ± 1 6 ± 9 MeV /c 2 , T ¥ ( 4360 ) = 94 ± 32 ± 1 3 MeV, and for the second one m Y (4660) = 
4669 ±21 ±3 MeV/c 2 , r r(4660) = 104 ± 48 ± 10 MeV. 
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